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Abstract. We have investigated the effect of the number of primordial binaries on the relationship between the total number 
of detected binaries within globular cluster and its collision rate. We have used simulated populations of binary stars in glob- 
ular clusters : primordial binaries and binaries formed through gravitational interactions. We show that the initial number of 
primordial binaries influences the relationship between the number of detected sources and the collision rate, which we find to 
be a power law. We also show that observing an incomplete sample provides the same results as those obtained with a com- 
plete sample. We use observations made by XMM-Newton and Chandra to constrain the formation mechanism of sources with 
X-ray luminosities larger than 10^' erg s"', and show that some of the cataclysmic variables within globular clusters should be 
primordial objects. We point out a possibly hidden population of neutron stars within high mass globular clusters with a low 
collision rate. 

Key words, binaries: general — globular clusters: general — X-rays: binaries 



^ 1. Introduction 



O 



X 



Globular clusters harbor many X-ray so urces (see for e xample 
the review of ROSAT observations by IVerbunt 11200 lb . Some 
of these sources were found to be foreground or background 
objec ts not related to the clusters (e.g. Verbunt & Johnston 
2000|). but others are c learlv associated with the clusters (e.g. 
Grindlav et al. Il2001a^. Using the latest X-ray observation fa- 



cilit ies, it has been shown through X-ray spe ctroscopy (see 
e.g. iRutledge et al. l2002HGendre et al. l2003bl). and bv optical 
follow-up observations (for example 'Grindlav et ai~ll2001allA 
[Pooley et al. 2002a; Webb et al. 2004) that the X-ray sources 
associated with the globular clusters are a mix of different bina- 
ries (cataclysmic variables, active bi naries, low mass X-ra y bi- 
nar ies) and millise cond pulsars (see I Verbunt & Lewin 1200 4). 

iHenon pointed out that a globular cluster made up 

of only single stars is unstable, and it will collapse in a few re- 
laxation times. The relaxation time is defined as the time when 
the exchange of energy between cluster stars due to close en- 
counters becomes significant for the cluster structure (see the 
review of Meylan & Heggie 1997). Because this time is short 
(a typical val ue for a galactic g lobular cluster is a few hundred 
million years, lDiorgovkill993l) compared to the globular clus- 
ter ages and because we observe many clusters that have not 
yet collapsed, several mechanisms have been suggested to ex- 
plain this : binary interactions, mass evaporation and black hole 
formation (Elson et al. 1987; Hut et al. ..1992) . As pointed out 



bv lElson et al. 1 lll987h . the mass evaporation should be most 
efficient at early times in the history of the cluster. Also, there 
is no firm evidence of a black hole in galactic globular clus- 
ters (Verbunt & Lewin 2004). On the other hand, simulations 
have shown that b inaries can prevent the c luster core from col- 
lapsing (see e.g. ^Baum gardt et al. ll2002h . leading to a stable 
configuration of the cluster, even if there is only a small sam- 
ple of binaries within the cluster This is due to the fact that 
binaries that interact with a cluster star can increase their bind- 
ing energy and provide the star with the excess kinetic energy 



ing energy ana proviae tne star witn tne excess Kinetic energy 
(lHeggiell975 ) . Thus , binaries can help to keep the cluster sta- 
ble fautetaLlI 19921) . and binary interactions within globular 
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clusters is a key point in studying the dynamical stability of 
globular clusters. Because most X-ray sources associated with 
globular clusters are binaries (or their offspring). X-ray obser- 
vations of globular clusters can provide strong constraints on 
binary formation. Here, we study these X-ray clues using the 
latest observations made by Chandra and XMM-Newton. 



This paper is organized as follows. In Sec.|2 we briefly dis- 
cuss the different ways a binary can form, and we explain the 
encounter rate within a cluster. In Sec.|3 we review the obser- 
vations made to date with Chandra and XMM-Newton as well 
as their main results. In Sec. |3 we study, from a theoretical 
point of view, what should be expected from the observations, 
using simulated samples. We compare the results from the sim- 
ulated samples to the observations in Sec.|5l 
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2. The encounter rate 

Binaries can form at the birth of a globular cluster (primor- 
dial binarie s), or later, by dynamical processes (dynamical bi- 
naries) (see IHuI et al. 1 19921 for a review about binaries within 
globular clusters). The number of dynamical binaries should 
scale with a parameter called the encounter rate (Clark 1975| 
rpabian et al. ll975UVerbunt & Hutll987tlVerbuntl20Q3l) . 
According to lVerburitI ( l2003l) . the encounter rate is : 

r = /.„cMv (1, 

In this equation, Ac,x represents the cross section of an in- 
teraction between species C and x, and n, represent the density 
of the species i. This value can be simplified with assumptions 
about the symmetry, equilibrium and gravitational potential of 
the cluster This gives the expression (see Verbunt 2003 and 
Hut & Verbunt 1983 for details): 

r-pi''-? (2) 

The encounter rate can thus be calculated using only the 
core radius value, r^., and the cluster central stellar density 
value po- If one of the hypotheses made in order to simplify 
this rate is false, then the encounter rate cannot be calculated 
with Eqn. |2 for this cluster co Cen provides a good example 
of such a limitation. Many sources in t his cluster are detecte d 
outside the inner part of the cluster jGendre et al. I l2003al) . 
This fact conflicts with the mass segregation: the most mas- 
sive sources (such as binari es) should gravitate to ward the cen- 
ter of globular clusters ( Me vlan & Heggie Il997l) . The equilib- 
rium state of a cluster implies that mass segregation has taken 
place (|Mevlan & Heggie 1997). Because mass segregation is 
not achieved in this cluster, the condition of equilibrium is not 
valid. 

These conditions imply use of this approximation only 
for globular clusters which are not collap sed (the cor e ra- 
dius is not defined in core-collapsed clusters, 'King '1962V For 
these collapsed clus ters, a more accurate method, used by e.g. 
iPoolev et al. 1 ll2003l) . is to compute numerically the integration 
of the Eqn. C] 

In order to arrive at Eqn. l2l lVerhiintl ll2003ll assumes: 

«c = Ccpo (3) 

In Eqn. |3j Cc is the compact object density fraction, com- 
pared to the cluster central density. This parameter Cc needs to 
be constant within the cluster (or at least within its central re- 
gion, where the integration of Eqn.^is carried out). However, 
comparing several encounter rates calculated using Eqn. |2l re- 
quires us to suppose that Cc has the same value in each cluster 
of the sample. If one calculates the encounter rate using Eqn.|2 
and assumes that F = pl^^r^, care should be taken when using 
this hypothesis. 

We will now focus on the X-ray observations of globular 
clusters, and compare the results obtained with several different 
encounter rate calculation methods. 



3. X-ray observations 

With the new observatories XMM-Newton and Chandra, the 
detection and identification of X-ray so urces related to glob- 
ular clusters is now common (see e.g. iGendre et al. Il2003at 
IPoolev et al. I l2002aV This has lead to the study of the bi- 
nary number versus the encounter rate by several author s 
jnendre et al. l2003hHPoolev et al. l2003tEke et al. l2003h . 
However, there is still the problem of the completeness of bi- 
nary samples (most of the faint binaries, such as active binaries 
have still not been detected in each cluster). Thus these au- 
thors studied either the number of quiescent low mass X-ray 
binaries with a neutron star primary within globular clusters or 
the whole sample of detected binaries, without any assumption 
about their nature. 

3.1. Low mass X-ray binaries within globular clusters 

Sources identified as a quiescent low mass X-ray binary with 
a neutron star primary (hereafter qLMXB'^'^) were classified 
as such objects due to th eir soft spectrurn and a luminosity of 
Ljc ~ 10^^ erg s ' (see e.g. lRutledge et al. l2002h . These objects 
with a soft spectrum define the sample used bv iGendre et al. I 
(2003b) and Poolev et al. (2003). We will now call them qNS. 
Due to their properties (soft sources with an X-ray luminos- 
ity around 10^^ erg s"'), qNSs are fairly easy to detect in each 
cluster observed by XMM-Newton or Chandra. These sources 
were thus thought to form a complete sample of objects. Using 
the XMM-Newton and Chandra observations, 'Gen dre et al. I 
(2003b), Poolev etal. (2003) and Heinke et al. (2003) claimed 
that qNSs within globular clusters were formed by stellar en- 
counters, and are thus not primordial binaries. 

Nevertheless, another kind of qLMXB'^'^ was detected by 
Wiinands et al. ( 2004). These objects have a power law spectra 
not compatible with th e soft spectrum of qNSs. According to 
IWiinands et al. I ll2004l) . these objects cannot be securely clas- 
sified using their X-ray color only because cataclysmic vari- 
ables have similar spectra. Thus, some sources classified as cat- 
aclysmic variables could be in fact qLMXBs JWiinands et al. I 
12004), and the qLMXB'*''^ sample should not be considered to 
be complete. 

3.2. Sources detected within globular clusters 

Using numerical integration of Eqn. [J IPoolev et al. I ll2003l) 
indicated that their whole sample of globular cluster X-ray 
sources was compatible with a power law relationship between 
the number of binaries and the encounter rate. Using the en- 
counter rate calculated with the simplified Eqn.|2l and another 
sample which meets the conditions required by using Eq n. |2l 
Gend^ (l2004l) confirmed this result. lHeinke et al1ll2003h used 
another method and fitted their sample with a law following 
the r^f/^^ relationship (a and /3 left free during the fit). Because 
the encounter rate is given by Eqn. ^ this method implies the 
same limitations as the ones indicated in deriving Eqn.|2l Their 
best fit values were respectively 1.0 and 1.3, and their sample 
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excluded the values of 1.5 and 2.0 (iHeinke et al. II2003I) . One 
should note that : 



(4) 



Thus, one can understand the 'Heinke et al. I (l2003l) work 
as a claim of a power law relationship between the number 
of X-ray binaries a nd the encounter rate given by Eqn. |2] 
iHeinke et al. I ( l2003l) have also studied the effect of the lumi- 
nosity limit variation or the spectral characteristics of their 
sample. They again found a correlation which can be inter- 
preted as a power law relationship. Moreover, they report a 
change in their best fitted values that can be interpreted as a 
change in the power law index. 

Three different groups, using different methods of collision 
rate calculation, have thus obtained the same result : the num- 
ber of detected sources located within globular clusters scales 
as a power law of the cluster encounter rate. Each group re- 
sults are independent ; ( iendre (2004) does not use the re- 
su lt and the enc o unter rate values fr om IPoolev et al.' J2003) 



"Heinke et al. ' ("2003"), and neither IPoolev et al. I ( l20'0 3) 
IHeinke et al. (2003) use the results from other groups. There 
is thus no problem when one compares these results as long 
as two different samples (from two different groups) are not 
merged without recalculating the encounter rate in a common 
manner. 

In this paper, we compare the relationship between the 
source number and the encounter rate. In Sec.l^we will use 
the encounter rate as a number, and we will not focus on how 
this number is obtained (supposed valid). In Sec. [S] we will use 
the sample of Gendre et al. (2003b) and thus his collision rate 
calculation method to discuss the number of qLMXB'*''^ within 
globular clusters. 

One should then explain why a linear relationship, observed 
for the qNS sample, appears to become a power law relation- 
ship for other sources. There are three possible explanations for 
such a behavior. First, the hypotheses made for the derivation 
of the collision rate are not valid for all the objects. Second, this 
could be due to the sample completeness. Finally, there could 
be a fraction of primordial binaries within the sample. We have 
examined the last two possibilities using simulated samples. 



(all luminosities have the same probability). These distributions 
were chosen to simulate two observational cases. The Gaussian 
distribution can deal with a rising luminosity function which is 
cut off at low luminosities by the instrument sensitivity. We 
have chosen the Gaussian parameters to be 4.0 x 10^' (mean) 
and 1.0x10^' (standard deviation) erg. s"'. The flat distribution 
deals with a unique population of object, where one should ex- 
pect that each luminosity has the same probability of appearing. 
The boundaries of this distribution have been chosen to be 10"^" 
and lO-'-' erg.s"'. 

To test the effect of completeness, we then excluded a part 
of the sample. Because we are dealing with an observational 
bias, we have excluded the faintest binaries of the sample. We 
have chosen to exclude 5, 10, 25, 50 or 75 % of the binary 
population. The flat and Gaussian distribution parameters are 
chosen arbitrarily, but because we exclude a percentage of the 
binary population, we can rescale the distribution parameters 
to any realistic parameters without changing the simulation re- 
sults. We then tried to find a power-law relationship between 
the collision rate and the number of remaining objects. The 
simulation is done 500000 times for both luminosity distribu- 
tions with a given threshold of binary exclusion (thus, a total of 
5 000 000 different simulations). At the end of each simulation, 
we compute the value of the power-law relationship index. At 
the end of all simulations, we compute the mean values of the 
indexes, which are listed in Tabled For both distributions, the 
power law index calculated is always compatible with unity. 

We have also looked for the effect of adding to the num- 
ber of dynamical binaries a random number value. This ran- 
dom number was chosen according to a Gaussian distribution. 
We then looked for a power-law relationship between the en- 
counter rate and this number of dynamical binaries added to a 
random value. The simulation is again repeated 500000 times 
for each of the Gaussian parameters listed in Table |2 The re- 
sults of these simulations are given in Table|2] in the 'popula- 
tion r column. 

The number of primordial binaries sho uld scale with the 
initial parameters of the cluster ( Murray et af Il99ll) . Because 
the encounter rate is based on current parameters, and because 
the evolution of the initial parameters depends also on param- 



4. Simulations 

We have found a linear relationship for a sample that we 
thought to be complete, but a power law relationship for an 
incomplete sample. The most natural explanation thus involves 
the completeness of the sample. To test this, we have used a 
simulation of a globular cluster sample. We have chosen a sam- 
ple of 33 clusters. Their encounter rates are fixed in order to be 
regularly spaced between two limiting values (0 and 100 using 
the conventions of Gendre et al. 2003b). First, we populate our 
sample with binaries using the linear relationship between the 
number of sources and the collision rate. We have chosen the 
number of binaries in a given cluster to be 0.5 x F. 

We then assigned a luminosity to each binary in each clus- 
ter We can use two diff'erent distributions: a Gaussian-like (i.e. 
there is a 'most probable' luminosity), and a flat distribution 



Table 1. Simulation of the effect of completeness. We have 
simulated two distributions of sources which obey the linear 
relationship given in Eqn. |2l We have then excluded a part of 
this sample and looked for a relationship. The power law in- 
dex is indicated in this table together with its error (errors are 
quoted at the Icr level). 



Sample 


Gaussian 


Flat 


excluded 


distribution 


distribution 


5% 


1.02 ± 0.03 


1.02 ± 0.03 


10% 


1.02 ± 0.05 


1.02 + 0.05 


25% 


1.02 + 0.06 


1.03+0.06 


50% 


1.00 + 0.08 


1.01 ±0.09 


75 % 


0.91 ±0.12 


0.92 + 0.12 
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eters external to the cluster (e.g. the galactic field), we should 
not expect the current encounter rate to be linked to the initial 
parameters. Because of this fact, they could be uncorrected. 
One can then assume them (and thus the number of primordial 
binaries) to be random numbers within the collision rate pa- 
rameter space. Thus, the random number added to the number 
of dynamical binaries could represent the primordial binariy 
number 

The assumed distribution of the primordial binariy num- 
ber (a Gaussian distribution) is somewhat arbitrary. Ideally, this 
distribution should take into account the evolution effects, the 
fact that not all primordial binaries are X-ray sources, and the 
fact that we do not detect all X-ray sources due to the lim- 
iting luminosity. Because most of these effects are now only 
partly understood, we have assumed that there is a most prob- 
able number of primordial binaries and modified the Gaussian 
standard deviation rather than testing other distributions (such 
as a flat or a cut-off' power law distribution). We have tried to es- 
timate the evolution effects, because gravitational in t eractions 
can d estroy primordial objects ( Pool ev et al. 
l200/ih . 

A primordial object is "destroyed" if it is disrupted or mod- 
ified by a dynamical encounter. Such an event occurs during a n 
exchange encounter (using the formulation of Verbunt 2003). 
According to .Verbunt ■ (.2003.) the exchange encounter rate is 
similar to the collision rate given in Eqn.^ The differences are 
that we now talk about the binary density (not the compact ob- 
ject density) and that the cross section depends on the binary 
size (not the stellar size). T he size of a bina ry is larger than that 
of a star (~3 stellar radii, Verbunt 2003); on the other hand, 
the number of primordial X-ray binaries (considered here, so 
~5 within each cluster) is less than the number of known com- 
pact obj ects within som e globular cluster (see e.g. the case of 
47 Tuc, iGrindlav et al. t2001a) . which is certainly less than the 
number of compact objects really present within these clusters 
(see Sec.|5}- Because we do not know precisely all of these pa- 
rameters, we have decided to use the same relationship to fix 
the number of primordial binaries dynamically destroyed and 
the number of dynamical binaries. We compared the number of 
primordial binaries which are present within the cluster (given 
by the random number generator) and the number of primordial 
binaries destroyed within the cluster If there are more primor- 
dial binaries destroyed than the number of primordial binaries 
present in the cluster, the total number of binaries is fixed to the 
value of dynamical binaries only. Otherwise, the total number 
of binaries is the sum of the number of remaining primordial 
binaries and the number of dynamical binaries. The search for 
a power-law relationship is finally repeated, and this simula- 
tion is again done 500000 times. Table Ogives the results of 
this simulation in the 'population 2' column. 

From Table |2l one can see that the power law index is no 
longer compatible with unity for either population 1 or popu- 
lation 2. The variation of the proportion of primordial binaries, 
simulated by varying the parameters of the Gaussian distribu- 
tion of the random number, always disagrees with a linear re- 
lationship. Moreover, as one can see in Table |2] the power law 
index decreases when the primordial binary fraction increases. 
These results are discussed in the next section. 



Table 2. Simulation of the effect of the presence of primordial 
binaries on a sample of dynamically formed binaries. We have 
simulated two distributions containing a mixture of sources that 
obey the linear relationship given in Eqn.|2land a random num- 
ber which simulates the primordial source number We have 
searched for a power law relationship within this distribution. 
Population 1 is a population where there is no dynamical de- 
struction, while population 2 is a population where primordial 
binaries can be destroyed by dynamical interactions. The power 
law index is indicated in this table together with its error (errors 
are quoted at the Icr level). The Gaussian parameters (in unit 
of sources) are the parameters of the distribution of the random 
number simulating the primordial binary number. 



Gaussian 


Gaussian 


Population 1 


Population 2 


mean 


sigma 






5 


3 


0.63 ±0.10 


0.72 ±0.10 


10 


3 


0.48 ± 0.06 


0.50 ± 0.06 


5 


5 


0.64 ±0.12 


0.70 ±0.13 


10 


5 


0.49 ± 0.09 


0.51 ±0.09 


20 


10 


0.36 ± 0.10 


0.29 ±0.10 



5. Discussion 

5.1. The qNS sample 

As already indicated, qNSs were thought to form a complete 



sample, and their number follow 


' a linear relationship with 


the encounter rate ( Gendre et al. 


2003b; 


Poolev et al. 


2003; 



then we also need to assume that the encounter rate calculated 
by Gendre et aL^ ( 2003b ) is valid, and thus that Cc is similar 
within each cluster. In other words, the neutron star density is 
constant within each cluster (the more massive a cluster, the 
more numerous the neutron stars are). 

Neutron stars can be detected as millisecond pulsars. As not 
all the millisecond pulsars within a globular cluster can be de- 
tected (due t o limiting sensitivity and the srnall ope ning of the 
cone ang le, iRansom lbOOSatlRansom et al. Il2003bl) . the num- 
ber of known neutron stars within one cluster should be taken 
as a lower limit to the real number of neutron stars within this 
cluster. The hypothesis that the neutron star density is constant 
cannot thus be verified. If this hypothesis is true, then a hid- 
den population of neutron stars is present, at least in massive 
clusters with no known neutron stars (such as M22). This pop- 
ulation could be composed of millisecond pulsars which have 
not yet been detected, but it could also consist of isolated cold 
neutron stars (which cannot emit X-ray light) or double neutron 
star binaries (which could be detected from their gravitational 
waves). 

We can also assume that Cr is not similar within each clus- 
ter, but then the relationship ^ Gendre et al. I ( 12003 b') found im- 
plies that this parameter cannot vary too widely, otherwise the 
collision rate might vary too widely from the computed values 
(and thus disrupt the relationship). We thus need to assume, 
again, a hidden population of neutron stars (despite the fact that 
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this population may be smaller than if the neutron star density 
did not vary at all). 

One may argue that the sample to consider would be the 
whole qLM XB^^ sample (qNSs, p lus the 'hard' qLMXB'^^ 
indicated by IWiinands et al. lE004l) rather than the qNS one. 
Because our current sample is then not complete, we can- 
not draw any conclusions. The 'hard' qLMXB^'^ EX01745- 
248 have very differe nt spectral properties than the qNS : 
IWiinands et aTl ( l2004 have shown that the canonical spec- 
tral model for so ft qNS (a neutron star atmosphere model, 

IPavlov et aTlllQQk.ZavUn et al 1996) cannot fit its spectrum. 

This may be due to a difference in the system properties, and we 
cannot rule out that 'hard' qLMXB'^'^ and qNS may be two dis- 
tinct kinds of qLMXB^"^, which could form in different ways. 
We also cannot rule out that these objects are the same kind of 
qLMXB, and that indeed our sample is not complete. But our 
simulations (see Tabled indicate that the non completeness of 
a sample does not change the correlation results or the value of 
the power law relationship (it only increases the error bars). 

5.2. The primordial binaries versus the dynamically 
formed binaries 

Our simulations have shown that there is no effect caused by 
the incompleteness of a sample (see Table [Q. Thus, any in- 
complete sample would allow us to study the correlation be- 
tween the number of objects and the encounter rate. Table |2l 
clearly shows that a population which is a mixture of primor- 
dial binaries and dynamically formed binaries has a correlation 
index which is less than one. The observations made by XMM- 
Newton and Chandra clearly indicate that the power law index 
of the correlation between the number of X-ray s ources within 
a cluste r and its collision rate is less tha n one. P oolev et al. I 
( l2002al) have found . 74 + 0.36, lGend're1 (>2004) found 0.5 + 
0.2 and lHeinke et al. I (l2003l) have shown that the value of one 
is excluded at the 90 % confidence level. This may indicate 
that the X-ray sources associated with the globular clusters are 
a mixture of primordial binaries and dynamically formed ob- 
jects. The variations of the exposant index es of the core radius 
and the central density, as reported bv Heinke et al. I (i2003) can 
be interpreted as a variation of the powerlaw index, and may 
then be due to a change in the luminosity limit or in the spec- 
tral properties of the sample which modifies the proportion of 
primordial binaries. 

At a limiting lumi nosity of L^^ = 4 x 10^" erg s"' (used 
by Poolev et al. 2003), the population is dominated by a mix- 
ture of qLMXB '^'^ and cataclysmic variables. If we assume (see 
above) that the 'soft' qLMXB'^'^ are formed by dynamical in- 
teractions, then the remaining objects ('hard' qLMXB'^'^ and 
CVs) should (at least for some of them) be primordial obj ects. 
One should note from the study of IHeinke et al. I ll2003h that 
their hard and luminous sample (their Fig. 5) is more compat- 
ible with a linear relationship (it is excluded at the 50% confi- 
dence level only) than their whole hard sample (their Fig. 6). 
This could indicate that the proportion of primordial binaries 
is less within the hard and luminous sample than in the hard 



whole sample. If so, this could be an indication of dynamical 
formation for some of these hard and luminous sources. 

6. Conclusions 

We have presented simulations of X-ray source populations 
and compare the simulated samples with the observations of 
XMM-Newton and Chandra in order to constrain the forma- 
tion mechanisms of binaries within globular clusters. Using the 
whole sample of detected X-ray sources within globular clus- 
ters, we have presented evidence to support the idea that glob- 
ular clusters contain primordial binaries, and that some cata- 
clysmic variables should be primordial objects. The growing 
sample of X-ray sources identified and thus associated with 
such kinds of objects will allow us to refine the X-ray samples 
and thus to improve these results. 
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